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Abstract 

Inelastic neutron scattering has been used to obtain a comprehensive descrip- 
tion of the absolute dynamical spin susceptibility x"(q) ^) of the underdoped 
superconducting cuprate YBa2Cu306.5 (Tc = 52K) over a wide range of en- 
ergies and temperatures (2 meV < fito < 120 meV and 5K < T < 200K). 
Spin excitations of two distinct symmetries (even and odd under exchange 
of two adjacent Cu02 layers) are observed which exhibit two different gap- 
like features (rather than a single "spin pseudogap"). The excitations show 

dispersive behavior at high energies. 

PACS numbers: 74.72.Bk, 61.12.Yp, 74.20.Mn 
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Typeset using REVT^ 



The magnetic excitation spectra of doped cuprates contain incisive and highly specific 
information about theories for high temperature superconductivity. Many models based 
on strong Coulomb correlations between charge carriers, for instance, predict that at high 
energies the spin dynamics should resemble those of the parent antiferromagnetic insula- 
tor. Until recently, neutron experiments on metallic and superconducting cuprates were 
confined to excitation energies below ~ 50 meV, larger than the superconducting energy 
gap but smaller than the intralayer nearest-neighbor superexchange J\\ ~ 100 meV which 
sets the energy scale for spin excitations in the undoped antiferromagnetic precursor com- 
pounds. Recent pioneering studies of the single-layer superconductor Lai.85Sro.i5Cu04 ex- 
tended these measurements to higher energies and established the presence of significant 
spectral weight at energies comparable to J\\ |l|J^. Analogous data for metallic and super- 
conducting YBa2Cu306+x (YBCO) have thus far not been reported. 

As YBCO is a bilayer system, such experiments can also answer questions about the 
nature and strength of the coupling between two directly adjacent Cu02 layers, a prob- 
lem of intense current research. As there are two Cu atoms per unit cell of YBCO, one 
generally expects the formation of bonding and antibonding electronic states. Conventional 
band theory predicts the formation of two Fermi surfaces in bands composed of bonding and 
antibonding states. However, theoretical as well as experimental [Q evidence indicates 
that these predictions do not generally hold for low-dimensional systems in the presence of 
strong correlations. Here we report neutron scattering measurements that evidence a new 
magnon dispersion-like behavior at energies of the order of J\\. Moreover, the spin excitation 
spectrum is characterized by two different gap-like features. These features provide funda- 
mental insights into the "spin pseudo-gap" phenomenon in metallic underdoped cuprates as 
well as essential new information about the interlayer interactions. 

Transitions between states of the same type (bonding-to-bonding or antibonding-to- 
antibonding) and those of opposite types are characterized by even or odd symmetry, respec- 
tively, under exchange of two adjacent Cu02 layers. In the cross section for inelastic magnetic 
neutron scattering, these transitions can be distinguished according to their distinct depen- 
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dences on the wavevector component L perpendicular to the Cu02 sheets The odd 

component of the spin susceptibihty displays a sin^(7rzcu-^) dependence whereas the even 
component has the complementary L dependence, cos'^{n zq^L). (Here, zq^ = 0.285 is the 
reduced distance between nearest-neighbor Cu spins within one bilayer, and the wavevector 
Q = {H, K, L) is measured in units of the reciprocal lattice vectors 27r/a ~ 2ti /h ~ 1.63A 
and 27r/c ~ 0.53A" ). In previous low energy neutron experiments in the metallic regime 
fSHTSi, only spin excitations of odd symmetry were found over the entire doping range. Does 
this reflect a fundamental selection rule, or is there a different (higher) energy scale associ- 
ated with the first type of transition? In the latter case, a comparison between the energies 
and susceptibilitites associated with the two excitation modes would yield quantitative infor- 
mation on the strength of the bilayer coupling, obviously an important parameter in models 
of high temperature superconductivity. This issue also bears directly on the interpretation 
of other experiments that are sensitive to this coupling. Our new measurements of high en- 
ergy spin excitations in underdoped metallic YBa2Cu306.5 have much better resolution and 
counting statistics than the initial studies on Lai.85Sro.i5Cu04 |]T],0]. Excitations of both odd 
and even symmetries are observed; the even excitations are characterized by a ~ 53 meV 



energy gap, the odd excitations show a gap-like feature at ~ 23 meV Both excitations 
exhibit large temperature dependences, and a dispersion-like behavior at high energies. 

Our sample was a high quality YBa2Cu306.5 single crystal of volume ~ 2.5 cm^ and 
superconducting transition temperature (midpoint) T^, = 52K. Its preparation and char- 
acterization have been described elsewhere |]13|; in particular, susceptibility measurements 



shown in Fig. 2 of Ref. |T3[ indicate a sharp superconducting transition (full width 5K) 
which rules out significant inhomogeneities in oxygen content. The neutron experiments 
were performed on the triple axis spectrometers INS (installed on a thermal beam) and 
INI (installed on the hot neutron source) at the Institut Laue-Langevin (ILL). The incident 
beam was monochromated by the (111) reflection of a flat copper crystal on INS, and by 
the (200) or (220) reflections of a vertically curved copper crystal on INI. On both spec- 
trometers, we used a graphite (002) analyzer with flxed vertical and adjustable horizontal 



curvatures. On INS, higher-order contamination was ehminated using a pyrolytic graphite 
(PG) filter on the scattered beam for different fixed final energies, Ef= 14.7, 30.5 and 35 
meV. On INI, we worked with a fixed final energy of Ef= 62.6 meV, and a nuclear resonance 
Er filter was placed on the scattered beam to suppress contaminations from higher energy 
neutrons. In all the measurements we scanned the wavevector Q while keeping the energy 
transfer constant. Two different scattering geometries were chosen in which wavevectors of 



the forms Q = {H,H,L) or Q = {3H, H, L), respectively, were accessible ||T2|. The results 
obtained on both instruments and in both geometries are in good agreement. 

Previous experiments [p|-[T5[| have established that the low energy magnetic cross section 
is peaked around the in-plane wavevector q2D = (fjf)? or H = K = ^ in reciprocal 
lattice units (r.l.u.). The cross section for odd spin excitations exhibits maxima at L^dd ~ 
1.7 + 3.5n (n = integer), and that for even excitations is maximum for Leven ~ 3.5n. We 
therefore performed constant-energy scans along Q = {H, H, Lodd) and {H, H, Leven), shown 
in Fig. 1. There are numerous optical phonon modes in the energy range covered by 
these data. However, scans in different diffraction zones as well as checks against phonon 



structure factor calculations |T2| established the magnetic origin of the peaks. The strong 
temperature dependence of the peak intensities (see below) is also inconsistent with phonon 
scattering. This methodology has previously been applied to neutron data at lower energies 
Il5|-|7|,|9|- |l2l , |l4| , and the results were found to be consistent with polarized beam experiments 



wherever such checks proved feasible PJT^. The temperature dependence of the uniform 
background presumably arises from multiphonon scattering. 

The scans were fitted to Gaussian profiles convoluted with the instrumental resolution 
function, corrected for the Cu magnetic form factor and converted to the dynamical 



spin susceptibility x" ^3 adjusting for the thermal population factor. Fig. 2 shows 
Xodd/even(^) = / '^q2DX" (q2D , ^odd/even , t^) / /c?q2D, the susccptibility averaged over the two- 
dimensional Brillouin zone in both odd and even channels. [We assume an isotropic q- 
dependence of the spin susceptibility around (7r,7r). This assumption is consistent with 
the data along (H,H) and a more limited data set along (3H,H).] While odd excitations 



are observed over the entire energy range probed by our experiments, an energy gap of 
^even ~ 53 meV exists for even excitations. Both Xodd('^) ^^w energies and the gap 
for Xeven('^) consistent with previous measurements in this doping regime |^,|^J^ which 
estabhshed a lower bound on the gap. Further, previous studies had given the neutron cross 
section in arbitrary units, but many quantitative models now require an absolute unit scale 
for x". We have therefore calibrated the magnetic intensity against the phonon spectrum, 
both against acoustic phonons at low energies and against an optical phonon at hu = 42.5 
meV according to a procedure discussed elsewhere [|12|. Both normalization procedures are 
in good agreement. 

The temperature dependences of both Xodd x'Len striking. At 200K, Xodd('^) 
shows a broad peak around 30 meV which sharpens and shifts to lower energy with decreasing 
temperature. The ~ 50 meV dip in Xoddl"^) lo"^ temperatures had not been observed 
before and may be related to the gap in Xevcn(^)- Much of the temperature evolution of 
Xoddl*^) takes place in the normal state; however, there is also an additional enhancement of 
the peak intensity at T^ which is presumably related to the magnetic resonance peak in 



the superconducting state dominating the magnetic spectra of more heavily doped samples 
(x > 0.9) |§-[TB|. While Xoddl'^) is little affected by temperature above fioj ~ 60 meV, at 
these energies Xeven('^) increases by almost a factor of two between 200K and 5K. In the 
same temperature interval Acven softens from ~ 59 meV to ~ 53 meV. This parallels the 
enhancement and softening of the peak in Xodd(^) lower energies. 

We now turn to the wavevector dependence of the spin susceptibility. At low energies 
(below ~30 meV), the scans are peaked around q2D with an intrinsic width of Aq2D ~ 0.2 
A~^ (FWHM). (Previous measurements with better q-resolution had indicated a flat-topped 
shape of these profiles [0). Fig. 1 shows that at high energies the peak broadens to 
Aq2D ~ 0.3 A~^ and disperses away from q2D = {\^\)^ so that a double peak structure 
emerges. For the sake of simplicity, we have fitted the q-scans to the sum of two displaced 
peaks only when a single peak did not give a satisfactory fit. The positions and intrinsic 
widths of the peaks in Xodd (^20, ^) and XevenlQao, ^) resulting from this procedure are shown 



in Fig. 3. 

It is instructive to compare these data to previous high energy measurements on the 
bilayer antiferromagnetic parent compound YBa2Cu306.2 and on the single-layer, optimally 
doped compound Lai.86Sro.i4Cu04. Clearly, the basic structure of the spin excitation spec- 
trum of YBa2Cu306.5 bears some resemblance to the spin wave spectrum of the YBa2Cu306.2 
i6| Jl7|j20| -p2| , with the ungapped odd (gapped even) excitations corresponding to acoustic 
(optical) spin waves, that is, in-phase (antiphase) precessions of localized spins in adjacent 
layers. The dynamical susceptibility of the acoustic spin wave mode is 

1 - 7^(q2D) 



where 7(q2D) = \ [cos^qxo) + cos(gy6)], and the quantum corrections for the spin wave 
velocity and the spin susceptibility are taken into account respectively as Z^. = 1.18 and 
Z^ = 0.51 following theoretical predictions p3[ in agreement with experiment ||12[. An 
analogous equation holds for even excitations, with cos^{tczcvlL) instead of the sin^ factor 
and a dispersion with a gap of 67 meV |^T|. The average of (1) over the 2D Brillouin zone is 



finite, 4:SZ^/ J\\fi%, and approximately independent of energy in the energy range probed by 
our experiment. For comparison this quantity is shown as the dashed line in Fig. 2 (as 10 
/i^/eV^ with J|| = 102 meV after quantum corrections [|l^). Clearly, the spectral weights of 
spin excitations in antiferromagnetic YBa2Cu306.2 and superconducting YBa2Cu306.5 are 
comparable. Finally, the spin wave dispersions of YBa2Cu306.2 are superimposed on the 
data of Fig. 3. Heuristically, the odd excitation spectrum of YBa2Cu306.5 can be fitted 
to a dispersive mode with a pseudo-gap of 23 meV |1^ and a dispersion of ~ 420 meVA, 
smaller than the antiferromagnetic spin wave velocity of 650 meVA [p!7tp2| . However, this 
analysis requires a very large damping parameter comparable to the gap. This pseudo- 
gap is characteristic of the normal state; at this doping level, the dynamical susceptibility 
is only weakly modified by superconductivity. Other features of the spin excitations of 
YBa2Cu306.5 such as the strong temperature dependence, the pronounced peak-dip structure 
at low temperatures and the energy- dependent broadening of the q-width are also very 



different from spin waves in YBa2Cu306.2- An explanation of these features presents a 
challenge to theories of the spin dynamics in the cuprates p^ . 

The reports on Lai.85Sro.i5Cu04 did not contain information about the temperature 
evolution of ^ind the measurements did not have sufficient resolution to establish the 

presence or absence of a dip feature in the spectrum. However, the general shape of 2D q- 
integrated x"i!^) iii deeply underdoped YBa2Cu306+x and optimally doped La2-xSrxCu04 
are substantially similar. Despite the different q-dependences of x"(q, cu) at low energies 
(four sharp incommensurate peaks in Lai.85Sro.i5Cu04, one broad commensurate peak in 
YBa2Cu306.5), both spectra are peaked near 25 meV. However, the maximum of x" P^r 
Cu02 layer is almost twice as large in absolute units in YBa2Cu306.5 (If given per formula 
unit, as in Fig. 2, the susceptibility of the bilayer compound exceeds the one of the single- 
layer compound by about a factor of 4). By contrast, fully oxygenated YBa2Cu307 has a 
much smaller normal state susceptibility and a sharp resonance peak in the superconducting 
state IHg]. 



In summary, we have reported the first observation of the even part of the dynamical 
spin susceptibility of the bilayer superconductor YBa2Cu306.5, an important part of the 
experimental description of spin fiuctuations in the cuprates. Both the odd and the even 
components are characterized by strong and unusual temperature dependences and by a 
dispersive behavior at high energies. Our observation of two different gap-like features in 
the dynamical spin susceptibility demonstrates that the spin-gap phenomenon, a hallmark 
of the underdoped cuprates, is more complex than previously thought. Our measurements 
of the spin correlations in absolute units over a wide range of energies, wavevectors and 
temperatures provide an excellent basis for theories of this phenomenon. 
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Figure Captions 

1. Constant energy scans of (a) Xodd(Q)^) through Q = (i, |,Lodd) and (b) XevenlQ'^) 
through Q = (|, ^,Leven), uieasured at 5 K (closed circles) and 200 K (open squares). 
Data have been obtained mostly on INS with a final energy of £'^=35 meV except at 
the energy transfer of huj =15 meV where Ef=14.7 meV. The scan for hu = 80 meV 
has been obtained on INI by rocking the sample around Q=(0.5,0.5,5.4) with £'^=62. 6 
meV. Data obtained on IN8 in different counting times have been rescaled to the the 
same monitor (monitor=2000, ~ 6 mn/point). Lines are the results of fits to a single 
Gaussian on top of a linear background, except at 80 meV where the line corresponds 
to two identical Gaussian lines displaced symmetrically on each side of if = 0.5. 

2. Odd (a) and even (b) spin susceptibilities at T=5 K, T=60 K (just above Tc) and 
T=200 K in absolute units (see text). Measurements using different final neutron 
energies Ef obtained on the two different spectrometers have been rescaled to the 
same units. The error bars do not include a ~ 30% normalization error. 

3. Spin excitation spectrum for odd (open symbols) and even (closed circles) excitations 
at 5 K. The open square indicates the energy of the maximum of the odd susceptibil- 
ity. The horizontal bars represent the full width at half maximum after a Gaussian 
deconvolution from the spectrometer resolution. The dotted lines correspond to the 
spin- wave dispersion relation in the insulating antiferromagnetic state with J[[ = 120 
meV and J±/J\\= 0.08 (without quantum corrections) |^ , ^ . 



8 



REFERENCES 



[1] K. Yamada et al, J. Phys. Soc. Jpn. 64, 2742 (1995). 
[2] S.M. Hayden et ai, Phys. Rev. Lett. 76, 1344 (1996). 

[3] See, e.g. D. Clarke, S.P. Strong and P.W. Anderson, Phys. Rev. Lett. 74, 4499 (1995). 
[4] See, e.g., H. Ding et a/., Phys. Rev. Lett. 74, 2784 (1995). 

[5] J.M. Tranquada, P.M. Gehring, G. Shirane, S. Shamoto and M. Sato, Phys. Rev. B 46, 

556 (1992). 

[6] J. Rossat-Mignot, L.P. Regnault, P. Bourges, P. Burlet, C. Vettier, and J.Y. Henry in 

Selected Topics in Superconductivity, Frontiers in Sohd State Sciences Vol. 1., Edited 
by L.C. Gupta and M.S. Multani, (World Scientific, Singapore, 1993), p 265. 

[7] J. Rossat-Mignot, L.P. Regnault, C. Vettier, P. Bourges, P. Burlet, J. Bossy, J.Y. Henry, 
and G. Lapertot, Physica C 185-189, 86 (1991). 

[8] H.A. Mook, M. Yethiraj, G. Aepph, T.E. Mason, and T. Armstrong, Phys. Rev. Lett. 
70, 3490 (1993). 

[9] L.P. Regnault, P. Bourges, P. Burlet, J.Y. Henry, J. Rossat-Mignod, Y. Sidis and C. 
Vettier, Physica C, 235-240,59,(1994); Physica B, 213&214, 48, (1995). 

[10] H.F. Pong, B. Keimer, P.W. Anderson, D. Reznik, F. Dogan, and I. A. Aksay, Phys. 
Rev. Lett. 75, 316 (1995). 

[11] P. Bourges, L.P. Regnault, Y. Sidis and C. Vettier, Phys. Rev. B 53, 876 (1996). 

[12] H.F. Fong, B. Keimer, D. Reznik, F. Dogan, and I. A. Aksay Phys. Rev. B. 54, 6708 
(1996). The units for x" used in this paper differ by a factor of 27r//|/3 from the ones 
used here. 

[13] H.F. Fong, B. Keimer, F. Dogan, and LA. Aksay, Phys. Rev. Lett. 78, 713 (1997). 



9 



[14] P. Bourges, L.P. Regnault, Y. Sidis, J. Bossy, P. Burlet, C. Vettier, J.Y. Henry, and M. 
Couach, Europhys. Lett, 38, 313 (1997). 

[15] P. Dai et ai, Phys. Rev. Lett. 77, 5425 (1996). 

[16] There is no unique way of specifying the pseudo-gap magnitude. Here, the spin pseu- 
dogap energy is identified as the energy at which dynamical susceptibihty is maximum. 
This definition then differs from previous conventions [Q. The normal-state spin pseu- 
dogap also differs from the true spin gap observed only in the superconducting state 
which occurs around 2-5 meV for x ~ 0.5 ||^,|^. 

[17] S. Shamoto, M. Sato, J.M. Tranquada, B. Sternlieb, and G. Shirane, Phys. Rev. B 48, 
13817 (1993). 

[18] x" refers to Tr(x"^)/3, where Xai3 is the imaginary part of the generalized susceptibility 
tensor. See S.W. Lovesey, Theory of Neutron Scattering from Condensed Matter, Vol 2, 
(Clarendon, Oxford, 1984), Chap. 8 p 39. 

[19] B. J. Sternlieb, J.M. Tranquada, G. Shirane, M. Sato, and S. Shamoto, Phys. Rev. B 
50, 12915 (1994). 

[20] J.M. Tranquada, G. Shirane, B. Keimer, M. Sato and S. Shamoto, Phys. Rev. B 40, 
4503 (1989). 

[21] D. Reznik, P. Bourges, H.F. Pong, L.P. Regnault, J. Bossy, C. Vettier, D.L. Milius, LA. 
Aksay, and B. Keimer, Phys. Rev. B 53, R14741 (1996). 

[22] S.M. Hayden, G. Aeppli, T.G. Perring, H.A. Mook, and F. Dogan, Phys. Rev. B 54, 
R6905 (1996). 

[23] J. Igarashi, Phys. Rev. B 46, 10763 (1992); A.J. Milhs and H. Monien, Phys. Rev. B 
54, 16172, (1996). 

[24] It is important to note that the similarities in the spectra of superconducting YBCO 



10 



with the antiferromagnetic undoped compound are not due to inhomogeneities in oxy- 
gen content. Apart from the sharp superconducting transition, The spectra of Fig. 2 
themselves are inconsistent with an energy independent admixture due to spin waves 
in undoped regions of the sample. Further, the low energy excitations, ~ 2 meV, are 
sharply suppressed below the superconducting Tq. (H.F. Fong et al, to be pubhshed). 



11 



FIGURES 



Odd 









2300 




2100 ■ 


o 


1900 1 






o 
E 


1700 • 


o 


650 - 


o 
o 


600 


04 


550 • 


CO 


500 


o 


450 ; 



>, 450 
1 350 ^ 
— 250 
1501- 
50 




0.2 0.4 0.6 0.8 1 
H (r.l.u.) 



12 



Even 



600 



550 - Ehl 

y - U 



500 ; 

450 • 



700 
650 
600 
550 
500 
450 
900 



5K ,T 
o 200K kf 




8001- 
700 
600 
500 



50meV 



0.2 0.4 0.6 0.8 1 
H (r.l.u.) 



13 




14 



1 50 



> 



E 



1 00 




C 
LU 



5 



